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SUMMARY 


This paper gives a brief consideration of rotary wing finite difference 
methods. The main concern is the specification of boundary conditions to prop- 
erly account for the effect of the wake on the blade. Examples are gxven of an 
approach where wake effects are introduced by specifying an equivalent angle of 
attack. An alternate approach is also given where discrete vortices are intro 
duced into the finite difference grid. The resulting computations of hovering 
and high advance-ratio cases compare well with experiment. Some consideration 
is also given to the modeling of low-to-moderate advance-ratio flows. 
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Hj, = tip Mach number 
R = blade radius 
r = radial station 

t = time 

Y = forward velocity 

= induced velocity 

X = x'/c, nondimensional streamwise coordinate 

Y = y’/Rs nondimensional spanwise coordinate 

Z = z*/R<S^^^, scaled nondimensional normal coordinate 

a = total angle of attack 
ttq = shaft tilting angle 

= angle of attack due to shaft tilting 
= angle of attack due to flapping 

Uj = angle of attack due to induced velocity 

3 = flapping angle 

Y = specific heat ratio 

6 = thickness ratio 

0^ = collective pitch angle 

0^ = twist angle 

y = advance ratio 
a = solidity 

= velocity potential 

= normal derivative of <{> 

= azimuthal angle 

to = angular speed 
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1 . Introduction 


The aerodynamic prediction tools available for routine rotor flow predic- 
tion and design have been limited to the simplest methods by the cost of compu- 
tation. However, current rapid decreases in the cost of computation are enabling 
a quiet, but fundamental reappraisal of our flow-prediction methodology. At 
present, integral methods (i.e., panel or lifting surface) with prescribed wakes 
are the most comprehensive methods being used in a routine fashion. However, 
these metliods are currently unable to treat nonlinear flow fields — most notably 
those associated with transonic flows. The best current means of treating such 
problems is with finite difference methods (FDM). 

The essence of the FDM approach is to use Taylor series expansions to 
approximate the partial differential equations by difference equations, which are 
then solved by implementing a variety of matrix methods. The method solves the 
flow region at a finite number of grid points. Proper flow resolution depends 
on placing a sufficient number of grid points at all significant flow regions. 
Unfortunately, a typical real rotor flow often contains so many important flow 
regions (the blade, the wake system, and the fuselage and acoustic propagation 
regions) that totally self-contained FDM solutions of complete rotor flows is 
not likely within the next five years. However, practical rotor loads, perfor- 
mance, and design work using FDM may be accomplished using isolated blade compu- 
tations, if we know how to properly specify our boundary conditions. The intent 
of this paper is to demonstrate various types of boundary condition specifications 
to perform FDM computations of various rotor flow conditions. For convenience, a 
small-disturbance potential equation will be used for this modeling. 

2. Equation of Motion and Boundary Conditions 

The simplest model equation which encompasses transonic flow effects is 
the classical transonic small-disturbance equation, which for a rotor has the 
form. 
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where F is a nonlinear flux term and is a function of (J)^ and 4)^. This equa- 
tion was derived under the assumption of small thickness and crossflow and high 
Mach number [1]. By the means of various mixed difference schemes, Eq. (1) has 
been shown to give good transonic solutions for unswept, nonlifting rotors for 
advance ratios as high as 0. 6 [2 ,3] . 

Equation (1) can be solved in three forms, 1) The steady form in which 
all time-derivative terms are ignored [4], In this case, high-speed hover solu- 
tions have been obtained using a classical relaxation scheme to solve the differ- 
ence equations. In ref. [4] the effect of the wake was obtained by including 
one vortex in the solution grid. This gave a solution which bore qualitative 
resemblance to a hover flow, but the model was too simple to be very accurate. 

2) The low-frequency solution in which all unsteady terms except are 

ignored [3]. This approximation is valid for flow variations which occur in 
several (say, 4) or more chord lengths of blade travel, which includes most 
rotor phenomena. Nonlifting computations using this model have been made using 
an alternating direction implicit (ADI) scheme and have compared very well with 
experimental data. However, good lifting solutions have not yet been obtained 
due to the lack of inclusion of a downwash model. 3) The fully unsteady case 
in which all unsteady terms are included. This case includes all high-frequency 
flows, the most notable of which is that resulting from the near passage of a 
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vortex. Perhaps the greatest hindrance to this problem is the lack of suitable 
experimental data. 

In general, Eq. (1) can be solved for simple flows; that is, the physics 
of compressibility no longer is a problem. However, the physical models required 
to accurately simulate real rotor flows (by FDM) have not been completely worked 
out or used. The rotor operates in the vicinity of a complex and often physi- 
cally close vortex system. Because we can only practically solve for the flow 
on an Isolated blade, the vortex system (as well as other flow features) must be 
accounted for by a careful specification of the boundary conditions. The 
simplest possible set of boundary conditions for an isolated lifting rotor must 
contain the following features (Fig. 1): 

1- The surface tangency condition . For the small disturbance equation 
this consists of equating the upper and lower surface slopes of the blade to the 
vertical flow component at a mean surface. 

2. The trailing vorticity sheet . This is an interior flow boundary 
condition — a sheet extending rearward from the blade, across which there is a 
specified potential discontinuity whose strength is the vorticity that has been 
convected from the trailing edge. This sheet is usually assumed to be undis- 
torted and coincident with a coordinate surface. 

3. The far-field boundary condition . The outer boundary surfaces are 
usually located sufficiently far from the rotor that free-stream conditions 

(}) = 0 or 4'n 0 can be specified there. It is important that the boundary 

condition used on an outer surface be compatible with any interior boundary that 
extends to this surface. For instance, compatibility with the trailing vorticity 
sheet requires that free-stream pressure be specified on the rear boundary. 

For a steady flow this corresponds to setting (})jj = 0. 

The inner radial surface boundary . To conserve mesh points, only the 
outer portion of the rotor flow is computed. Since the inner rotor flow is two 
dimensional away from the tip, it suffices to use strip theory (i.e., ((>nn = 0) 
on the inboard computational plane. 

Although the above boundary conditions are all that are required for fixed 
wing work, they are incomplete for the rotor situation because the downwash 
induced by the wake system is not accounted for (except for that flow induced by 
the vortex sheet included in the grid) . The simplest possible way to include 
this downwash is to modify the surface tangency condition by prescribing an 
effective angle of attack. An example of this approach is given in the follow- 
ing section. 

3- Computation of High-Speed Advancing Blade Loads. Comparison with Experiment . 

An ideal situation in which to employ the angle-of-attack approach is for 
high advance-ratio flight (say, greater than 0.3) where the induced downwash is 
small. Furthermore, there exists a good body of model blade surface pressure 
data for this case [2,3,5]. This data (obtained at ONERA using 2 and 3 bladed 
rotors with instrumented removable tips) covers a number of high advance-ratio 
conditions ranging from nonlifting to moderately loaded (C^/a = 0.06) cases. 

All previous computations of these cases have been only for the nonlifting cases 
using the low-frequency form of Eq. (1). To perform some first computations of 
the lifting cases, the Induced incidence has been calculated using the Drees 
downwash model [6]. This simple model neglects all discrete blade effects (such 


4 



3 .S tli 6 tip vorticBs) and computes downwasli only at tlie irotoir center and 3/4 
radius (the downwash elsewhere is taken as a linear interpolation or extrapola- 
tion from these two points). The simplifications of this model are such that it 
was not deemed necessary to modify this downwash by the vortex— sheet— induced 
velocities when it was used in the finite difference code. 

The unsteady calculations have been performed on an isolated blade with 
an angle of attack prescribed along the blade for each azimuthal location by. 

a(r/R,iJ^) = 0^ + 0^(r/R - 0.7) +01 + 0^+03 


with 


= tan ^[V sin cos + V cos Uq sin y\>)] 

= tan"^{[-rf2(d3/di|;) - V cos cos ^ sin 3] /(fir + V cos sin \|^) ] } 
= tan“^[-V^ cos &/(Qr + V cos aq sin ij;) ] 
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The experimental values are used for the rotor shaft only (otg) , the col- 
lective pitch angle (0^), the blade twist angle (0^) » the flapping angle (3). 

An example showing the evolution and relative magnitude of the different terms: 

0c + 0^(r/R - 0.7); a^; 02 ; and is given in Fig. 2. It is seen in this 
figure that although the induced angle is the smallest single contributer to the 
inflow, it is nevertheless significant in the first quadrant where the other 
terms nearly cancel each other. 

Two experimental cases have been computed corresponding to the following 
conditions; 
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First case: 


C^/a = 0.05, V = 76 m/sec, fiR = 210 m/s 


Second case: 


C^/o = 0.0665, V = 81.4 m/s, OR = 210 m/s 

^ The theoretical pressure distributions for the azimuths \li = 30°, 60°, 90°, 
120®, 150®, and 180® are compared with the experimental results in Figs. 3-5 
(Cip/cT — 0.05) and Figs. 6—8 (Cj/o = 0.0665). The evolution of the pressure 
distribution is fairly well predicted, including the development of a shock wave 
for \p = 90° and 120®. For the second case, the computed shock wave tends to be 
stronger and located further aft of the experimentally obtained shock. 

In general, the discrepancies between experiment and calculation tend to 
be greatest in the first quadrant. These differences between experiment and 
calculation are more easily seen if we display the time history of the pressure 
coefficient at several points on the blade surface. Such time histories are 
shown for the spanwise section, r/R =0.9, for Cj/o =0.5, and for 0.665 
(Figs. 9 and 10, respectively). The differences seen in the first quadrant are 
almost certainly a result of inaccuracies in the very simple Inflow model. 

(Recall that the induced velocity is much greater in the first quadrant than in 
the second.) The most significant error in the second quadrant is in the shock 
location. This may also be a residual effect of the induced velocity. However, 
an equally likely possibility is that we are seeing the effects of neglecting 
the boundary layer. It has been shown, by a number of calculations with boundary 
layer models included, that the effect of viscosity is to weaken the shock and 
move it upstream [7]. 

Overall, the results obtained by these first lifting calculations of the 
flow around a helicopter blade in forward flight (y > 0.3) are in good agreement 
''^ith the experimental results. However, the use of a more sophisticated model 
for the determination of the induced velocity is certainly necessary for the 
lower advance ratios . 

Lifting-Rotor Computations by Intra-Grid Vortex Insertion. Steady Cases . 

Transonic computations using an angle— of— attack condition cannot .be valid 
when the vortex is too close to the blade surface. This is because the flow 
nonlinearity demands that the effect of the vortex be computed at all points in 
the flow, and not just at the blade surface. To handle this situation, one can 
conceive of solving for the perturbation potential about the vortex— induced 
velocity field. However, a simpler approach is to modify the existing vortex 
sheet logic referred to in Section 2. That is, a vortex line (of constant 
®brength for the hover case) can be inserted into the grid by specifying a sur- 
face across which there is a potential discontinuity whose magnitude equals the 
*^^®Lred vortex strength. The edge of this sheet is the desired vortex line. 

For a constant-strength vortex the sheet is merely a branch cut and its orienta- 
tion should have no effect on the solution. Such a vortex is shown in Fig. 1, 
represented by a vertical sheet. 

Hover flight is a good case to test the approach described previously 
because the vortex strength is constant. Furthermore, the vortex geometry is 
^^Lrly simple and can be reasonably represented by straight lines parallel to 
the flow in the near-blade region. For present purposes, a number of steady 
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finite difference computations have been made of hover flows using the relaxation 
code of Ref. [4]. In the current computations, two vortices (represented by 
vertically oriented sheets) are located beneath the blade. However, this does 
not yet account for the total effect of the entire wake system. We assumed 
therefore that the more distant elements of the wake system can be accounted for 
by a modified angle of attack. The modification consists of computing the angle 
of attack (using the Biot-Savart law) only from those elements of the wake system 
not contained in the finite difference grid. Figure 11 shows the computed load 
distributions on hypothetical twisted and untwisted methods. The FDM results are 
compared with those of the lifting-surface code, HOVER, described in Ref. [8]. 

Both codes use the identical far-wake model [8] and assume the vortex strength to 
equal the maximum blade-bound vorticity. The comparison of the two methods is 
quite close. The differences that do exist could stem from a number of causes 
including vortex curvature and boundary conditions. 

There also exists a convenient body of experimental data [9] by which to 
test the previously mentioned wake treatment. This data consists of extensive 
surface pressure distributions and vortex geometry measurements for a high-speed 
hovering rotor (untwisted blade, NACA 0012, >R=6.0). The vortex and blade 
geometry is the only required input for computing the rotor flow. A comparison 
of the experimental lift distribution with the finite difference and lifting- 
surface methods is shown in Fig. 12. At this rotor speed (M<j» = 0.6), the 
lifting-surface and finite difference results compare closely. However, there is 
a greater difference near the tip than occurs in the previous high aspect-ratio 
results. This is almost certainly a result of there being no wake curvature in 
the finite difference grid. Both the finite difference and lifting-surface 
results overestimate the maximum lift. Overall, however, the comparison with the 
data is good. Also included in Fig. 12 is a finite difference computation using 
only an angle-of-attack boundary condition. For this computation the angle of 
attack was obtained from momentum theory. It is not surprising to find that the 
resulting lift is seriously in error at the tip, due to the lack of discrete blade 
effects in the momentum model. However, the inboard results are good for this 
case. This indicates that such a simple downwash model may be generally useful 
for specifying the inboard flow regions. 

The comparison of lifting surface, FDM results, and experiment is espe- 
cially interesting at high tip speeds. Figure 13 shows such a comparison at a 
tip speed of Mj = 0.877. Since the lifting-surface method is linear (as are 
almost all integral methods) it is not surprising that this method greatly over- 
estimates the lift in the tip region. However, the FDM computations compare 
well with the experimental lift distribution. Figure 13 also shows a comparison 
of the measured pressure distributions with FDM results. The comparisons are 
favorable even to the extent of properly predicting the correct shock locations. 
This is a little surprising when one considers that the computed pressure per- 
turbations exceed the assumed limits for small perturbation theory. Indeed, the 
inaccuracies seen in the leading-edge region undoubtedly stem from small distur- 
bance errors. In addition, no boundary layer corrections have been made. 

All operational helicopters operate at low enough tip speeds that linear 
integral methods will always suffice for hover work. However, in demonstrating 
the ability to easily duplicate integral results by FDM and to effectively apply 
these to high speed, we will have performed an essential validation in prepara- 
tion for general forward-flight computations. 
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5. Vortex-Induced Effect in Forward Flight . 


The extension of the previous mode of vortex representation to unsteady 
cases has geometrical, numerical, and practical problems associated with it. To 
begin, the wake for these cases can no longer be represented by a straight, 
unmoving vortex of constant strength. The sheet defining the vortex system is 
now a physically real vorticity sheet (for a tip vortex of varying strength) and 
cannot be given the arbitrary treatment of the hover cases. The tip vortex 
assumes a wide range of orientations with respect to the blade and can be nearly 
parallel to it in certain encounters. Clearly, a general forward-f light flow 
treatment requires the fully unsteady form of Eq. (1). Although there are 
methods to solve the fully unsteady Eq. (1), there have not yet been many experi- 
ments performed to verify these methods. A considerable problem is that there 
currently exists little usable wake and loading data -by which to evaluate new 
forward-flight methods at the lower advance ratios (say, below 0.3). Therefore, 
at present the best approach for computation work is to calculate model problems 
which contain the qualitative effects of blade/vortex interactions in the simplest 
possible way. 

One such model is the two-dimensional interaction of an airfoil with a 
moving vortex. This model is particularly easy to implement because the vortex 
can (as in the hover case) be represented by an arbitrarily oriented sheet 
across which a potential discontinuity is specified (Fig. 14) . The edge of the 
sheet (the vortex) is currently stepped through the grid in a straight line at 
the speed of the undisturbed free-stream velocity. The strength of the vortex 
is given as an effective lift coefficient, of an airfoil having the same 

circulation as the vortex. The technique for solving the fully unsteady Eq. (1) 
in its two-dimensional form is a simple ADI scheme given in Ref. [10]- The 
effect of unsteadiness in these computations is shown in Fig. 15, which compares 
blade surface pressure distributions for a fixed and moving vortex. These dis- 
tributions are shown when the vortex is at the leading-edge, mid-chord, and 
trailing-edge chordwise locations and at 0.96 chords below the blade (Clv = 0.1). 
Although the steady case clearly shows the effect of the vorticity there is no 
apparent disturbance for the unsteady case. This result undoubtedly reflects 
the fact that the vortex exerts no net force on the fluid in the unsteady case. 

In fact, in order to get a sizable unsteady effect on the surface pressures when 
the vortex is beneath the profile, it is necessary to bring the vortex much 
closer to the blade and increase its strength. Figure 16 shows the lift varia- 
tion for a blade with a vortex of strength C^y =0.4 whose path is 0.26 chords 
beneath the blade. As the vortex passes beneath the blade it is seen to signifi- 
cantly enlarge and strengthen the supercritical flow region on the bottom surface. 
However, the effect quickly disappears after the vortex has passed the shock. It 
should be recognized that this latter computation could contain a serious error 
because it does not currently allow the vortex to follow the flow (which varies 
greatly from the free stream close to the airfoil). This error should cause an 
exaggeration of the effect of the vortex on the airfoil loads. It seems likely 
therefore that these close high-frequency interactions will have a far smaller 
effect on the loads than quasi-steady theory would indicate. This does not say, 
however, that the effect will be unimportant. The determination of the magnitude 
of this blade-vortex interaction must await a good experiment and more detailed 
theoretical modeling. 

5. Concluding Remarks 

The increasing availability of good FDM rotor codes and cheaper computers 
will ultimately enable the routine simulation of the nonlinear unsteady 
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three-dimensional flow about rotors. Initially, these simulations will be iso- 
lated blade computations whose specified boundary conditions account for the 
effect of the wake and interference effects. This paper demonstrates two 
approaches to the treatment of boundary conditions which result in realistic FDM 
computations of both hover and high advance-ratio (p > 0.3) rotor flows. For the 
latter case (computed with a low-frequency unsteady code) discrete vortex 
effects are specified as an angle-of-attack variation. The inflow model is very 
simple and can be improved, although good comparisons with data have been 
obtained. Even better comparisons have been made for high-speed hover computa- 
tions (using a steady code). However, for this case it is necessary to use a 
more exact boundary condition by actually inserting some of the tip vortices 
into the finite difference grid. In addition, a modified angle of attack (com- 
puted from those wake elements not included in the grid) is also required. For 
the low-to-intermediate range of advance ratios, the problem is more difficult 
because the wake has a complicated structure which cannot be ignored. Neverthe- 
less, demonstration computations (using a fully unsteady code) have been made on 
a two-dimensional blade/vortex interaction problem which contains the basic 
numerical (if not all the geometrical) features of the problem. 

The numerical techniques used in this paper are not new but they have not 
yet seen the full realization of their capabilities for rotor loads, performance 
prediction, and design work. The hover and high simulations performed in this 
paper may be the first truly realistic FDM rotor simulations. There seems to be 
little technical reason for reticence in the use of FDM techniques and we trust 
that the economic reasons will rapidly diminish. 
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Boundary conditions for the rotor potential flow problem. 






Cjh = 0.05, Vq = 76 m/sec, wR = 270 m/sec 



Fig. 2. Total angle of incidence calculated by Drees downwash model. 
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Fig. 5. Comparison of measured and computed chordwlse pressure distribution at 
different azimuthal angles, Cj/o = 0.05, r/R = 0.95. 
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Fig. 7. Comparison of measured and computed chordwise pressure distribution at 
different azimuthal angles, C^/a = 0.0665, r/R = 0.90. 
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Cj/a = 0.0665, Vq = 81.4 m/sec, , coR = 210 m/sec 
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,fi = *?n° 


r/R = 0.95 


1.0 0 


Fig. 8. Comparison of measured and computed chordwise pressure distribution at 
different azimuthal angles, Cj/a = 0.0665, r/R = 0.95. 





1.5 


Cj/a - 0.0665 , Vq - 81.4 m/sec, coR = 210 m/sec 




Fig. 10. Comparison of measured and computed pressure coefficient at different 
chordwise stations, Cr^lo = 0.0665, r/R = 0.9. 
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AR = 6.0, Mt = 0.6. 0g = 8° 


FINITE DIFFERENCE (ROTRAN) 

PANEL METHOD (HOVER) 




Fig. 11. 


Comparison of computed 


loadings by panel and 


finite difference method. 
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UPPER SURFACE » FINITE ° UPPER SURFACE \ 

LOWER surface/ DIFFERENCE • LOWER SURFACE 



Fig. 12. Comparison of measured and computed chordwise pressure distribution at 
Mrj, = 0.877, r/R = 0.89, and r/R = 0.96, respectively. 


22 




Mj = 0.877, AR = 6.0, 6 ^ = 8 ° 

UPPER SURFACE ) FINITE O UPPER SURFACE ) 

LOWER SURFACE ) DIFFERENCE • LOWER SURFACE ) 




Fig. 13. Comparison of measured and computed chordwise pressure distribution at 
Hj, = 0.877, r/R = 0.5, and r/R - 0.68, respectively. 
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SMALL DISTURBANCE BOUNDARY CONDITIONS 
FOR 2-D WING/VORTEX PROBLEM 



Fig. 14. Boundary conditions for the blade/vortex interaction problem. 
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problem. 
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Fig. 16. Unsteady flow computations for the blade/vortex interaction problem. 
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